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FOREWORD 


The work described in this report was performed in the Fire Control Formulation Branch 
(K41), Submarine Launched Ballistic Missile (SLBM) Research and Analysis Division, Strategic 
and Space Systems Department, and was authorized under Strategic Systems Program Office Task 
Assignment 36401. This work was necessitated by the need to formulate an exact method of 
computing geodetic latitude and altitude whose chief attribute is the absence of singularities in the 
results. 


This report is a refined version of NSWC TR 85-85, dated April 1985 and revised May 
1986. NSWC TR 85-85 should be discarded. 

This report has been reviewed and approved by Davis Owen; Johnny Boyles; Carol Rose, 
Head, Fire Control Formulation Branch; and Dr. David Lando, Head, SLBM Research and 
Analysis Division. 


Approved by: 



R. L. SCHMIDT, Head 

Strategic and Space Systems Department 
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INTRODUCTION 


A simple and efficient model is presented for calculating the exact geodetic latitude and altitude 
of an arbitrary point in space, given the coordinates of that point. The model assumes the Earth to 
be an oblate spheroid; i.e., an ellipsoid of revolution whose semimajor axis ‘a’ is the radius of the 
circle described by the equatorial plane, and whose semiminor axis ‘b’ is a line joining its center 
and one of its poles. The point in question can be either inside or outside the Earth ellipsoid, 
excluding a small open region bounded by a prolate spheroid that is concentric with and contained 
within the Earth ellipsoid. The cardinal merit of this model is the absence of singularities in the 
results. 


FORMULATION 


Choose an Earth-fixed Cartesian coordinate system whose origin coincides with the center of 
the Earth ellipsoid. The unit vectors 1 , j, and fc coincide with the x, y, and z axes, respectively. 
The +z axis points in the direction of the North Pole. The +x axis is the line of intersection of the 
equatorial plane with the plane of zero longitude. The +y axis completes a right-handed coordinate 
system. 

An equation of the ellipsoid in this frame is 

^rr+—? = 1, where r = Jx 2 + y 2 . ^ 

a 2 b 2 

Let P(xo, yo, zo) be the coordinates of the given point. It is desired to find the nearest point 
P'(x,y,z) to P on the surface of the ellipsoid. 

The slope of a normal to the ellipsoid at any point on its surface is given by 

1 a 2 z r 2 z 2 

--j— = —r-, where —r + —r = 1. (2) 

dz b 2 r a 2 b 2 w 

dr 

Therefore, the slope of a normal from P is 

z 0 z a z , / 2 j 2 " 

— -= - 5 -, where r 0 = ^x 0 + y 0 ; 

r 0 -r b 2 r 

i.e., 

(r 0 -r)a 2 z = (z 0 -z)b 2 r, 
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or 


r 0 a 2 z = r{(a 2 -b 2 )z + b 2 z 0 }. 

Squaring both sides and expressing r in terms of z, 

a 2 b\V=(b 2 -z 2 ){(a 2 -b 2 )z + b 2 z 0 } 2 . 

Writing the above equation in descending powers of z, the following quartic in z is obtained: 

(a 2 - b 2 )V + 2b 2 (a 2 - b 2 )z 0 z 2 + b 2 {a 2 r 0 2 + b 2 zj - (a 2 - b 2 ) 2 }* 2 - 2b‘(a 2 - b 2 )z„z- b‘zj = 0 . 

Since z is either positive or negative and sign z = sign Zq, the above quartic can be expressed in 
terms of z . The result is 


z| 4 + 2p|z| 3 + q|z| 2 - 2pb 2 |z| - p 2 b 2 = 0 , 


(3) 


where 


b 2 z, 


0 _ rol 


P a 2_ b 2 e /2 ’ 


(4) 


b 2 {a\ 2 + b 2 4-(a 2 -b 2 ) 2 } 

q= - V ^) 2 — 


_2 .2 . 

p - b + 


2/2 ’ 
e e 


(5) 


where 


e 2 = l-^-,e' 2 =4-l 

a b 2 


Two of three powerful theorems in the Theory of Equations will now be employed to expose 
the nature of the roots of Equation (3). The third theorem will be needed later on. These theorems 
are found in the reference 1 . 


I. An equation f(x) = 0 cannot have more positive roots than there are changes of sign in f(x), 
and cannot have more negative roots than there are changes of sign in f(-x). 

n. Every equation that is of an even degree and has its last term negative has at least two real 
roots, one positive and one negative. 

III. Every equation of an odd degree has at least one real root whose sign is opposite to that of 
its last term. 


By observing Equations (3), (4), and (5), the following conclusions are deduced: 


1 Hall, H.S. and Knight, S.R., Higher Algebra . Fourth Edition, MacMillan and Co., 
London, 1927. 
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(i) Since the last term of Equation (3) is negative, there are at least two real roots of opposite 
sign (using II). 

(ii) Since there is only one change of sign in Equation (3), there is at most one positive root 
(using I). 

From (i) and (ii), it is seen immediately that Equation (3) has exactly one positive root, and this 
is the root that is sought. 

The solution of Equation (3) is effected by a standard method known as Ferrari’s method 1 . In 
the equation 


|z| 4 + 2p|z| 3 + q|zf - 2pb 2 |z| - p 2 b 2 = 0 , 

add to each side (clzl+d)2, the quantities c and d being determined so as to make the left-hand side 
a perfect square; then 

|z| 4 + 2p|z| 3 + (q + c 2 )|z| 2 + 2(cd - pb 2 )|z| + d 2 - p 2 b 2 = (c|z| + d) 2 . 


Suppose that the left side of the equation is equal to (Izl 2 + plzl +1) 2 , then comparing the 
coefficients, 


p 2 + 2t = q + c 2 , pt = cd - pb 2 , t 2 = d 2 - p 2 b 2 . 
Eliminating c and d from these equations, 

p ! (t + b ! ) 2 = (2t + p 2 - q)(t 2 + p 2 b 2 ), 


( 6 ) 


or 


or 


2t 3 - qt 2 - p 2 b 2 (q - p 2 + b 2 ) = 0 , 


t) 2 b 2 r 2 

2t 3 - qt 2 - P r ° = 0, using Equation (5) 
e e 


or 


2 2 2 


2t 3 - qt 2 - —= 0 


(7) 


Applying Theorems I and EH to Equation (7), we see that Equation (7) has exactly one positive 
root. The solution of Equation (7) is accomplished by a standard method known as Cardan’s 
Solution 1 . Eliminating the t 2 term in Equation (7) by malting the substitution 
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f nV ( nV aW 

2C+ l - q,,+ § 

\ 6J \ 6J e 


where 


t' 3 +Ft' + H = 0 , 


a 2 n 3 a 2 rV 

p__ t. pp__3_ _ 

12’ ~ 108 2e' 8 


The solution of Equation (8) is 


, H 2 F 3 H 
1 “hi + 27 2 


H 2 F 3 H 
4 + 27 + 2 


Equation (9) is valid provided that 


We have 


where 


£l + ^>o. 

4 27 


n 3 aW 9 

H = ——-i^ = -4(q 3 +2P), 

108 2e 8 6 3 V ' 


108 2e' 


r . 27a r 0 z 2 
e' 8 


H 2 F 3 1 


T + f7^( q!+2p ) 2 -i^4 p ( p+q! ) 


»l + ^>o, 

4 27 


P + q 3 > 0 , 

which is satisfied when q > 0; i.e., when a 2 r 2 + b 2 z 2 ^ (a 2 - b 2 ) , which represents all space 
excluding the open region bounded by the ellipsoid 
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2 2 
-&_ + _5l = i 

aV b 2 e' 4 


Since ae 2 < be' 2 < 43 km, this region is of no practical interest. Therefore, q ^ 0 will be the 
constraint. Applying Theorem I to Equation (7) with -t substituted for t, we see that Equation (7) 
has no negative roots. Hence, Equation (7) has only one real root and it is positive. 


H*_ 

I 4 + 27 2 " 6 3 


where 


Similarly, 


where 


/p(p+q ! )+^(2P+q ! ) = ^-{2P+q 3 + 2 1 / p ( p +q 3 )} 

=4(i/p+q T +V p ) 2 . 


Ih 2 

F 3 


[i 

+ 27 

2 J 


Q=(V p + q 5 +Vp)’ • 


\W_ + P + H 3 
4 + 27 + 2 , 6 


Q' = (VP + q 3 -Vp) 3 . 


(Note that Q' = —. However, putting Q' in this form is undesirable as it introduces a singularity 
when Q = 0, which occurs when either both q and r 0 equal zero or both q and Zq equal zero.) 
Substituting these values into Equation (9), 

t'=7(Q+Q')=><=7(q+Q+Q')- 

O O 

It will now be shown that q < 2t < q + b 2 , where the left-hand equality occurs when either 
r 0 = 0 or Zo = 0 and the right-hand equality occurs when ar 0 = blzj,!. 


Rewriting Equation (7) in the form 


2 22 
3. r n z n 


t 2 (2t-q) = ^, 
e 


the lower bound of t is deduced at once. Now let 
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where t > q. Then 


3 2 4 ar 0 Zp 

g = x 3 -qt 2 -^ 

e 


= x(3x - 2q) > 0 . 


Thus, g is a monotonically increasing function of T for all x > q. Since 2t > q and g = 0 when 
x = 2t, it follows that g > 0 when x > 2t. Putting x = q + b 2 in the expression for g, 


i=, rt -=b ! (q + b=) 2 -^ = b 2 (p^-^l- 


2„2 1 

4b p r 0 _ . 

— r u , using Equation (5) 
e e 


=b 2 p 2 --^ *o 

v ee J 


r 0 r o Z o 

2 t <q +b 2 , equality occurring when p =— 7 ; i.e., when— = ——, or when ar 0 = b|z 0 |, in 

cc c c 


which case, 


q = 2 p 2 -b 2 , t = |(q + b 2 ) = p 2 . 


Solving for c and d from the system of equations ( 6 ), 

c = -^p 2 -q + 2 t, d = -(t + b 2 ) = Vt^+pV . 
c 

Since 2t > q, it follows that c > p. Now, (lz|2 + plzl + 1) 2 = (clzl + d ) 2 implies 
Izl 2 + plzl +1 = ± (clzl + d), from which the following two quadratics in Izl are obtained: 

|z | 2 + (p-c)|z| + t-d = 0 , 


|z| 2 +(p + c)|z| + t + d = 0 . ( 12 ) 

It has already been shown that Equation (3) has exactly one positive root. Furthermore, since 
t 2 — d 2 = -p 2 b 2 < 0 and d > 0, it follows that t-d<0<t + d. Now, t - d is the last term of 
Equation (11) and has just been shown to be negative. Hence, by applying Theorem II to Equation 
(11), we see immediately that Equation (11) must contain the required positive root. Hence, 
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-(p-c) + V(p-c) 2 -4(t-d) 
2 


= ~2 v ~ P + v 2 P 2 “ 9 “ 2t ~ 2pc + 4-i/t 2 + p 2 b 2 J . 

This solution is free of singularities. It is easy to show that lim|z| = 0 and lim|z| = b . We have 

z 0 —>0 r 0 —>0 

lim p = 0, lim P = lim P = 0 => lim Q = lim Q' = lim Q = lim Q' = q 

z 0 -»O z 0 -»0 r 0 —>0 z 0 -»0 z 0 —>0 r 0 -»0 r o -»0 


=> lim t = lim t = — => lim c = 0, lim c = p. 

z 0 -*0 r 0 —>0 2 z q-> 0 ig —*0 




To compute the geodetic latitude and altitude at P'(x,y,z), it is desirable to introduce a 
geometric term, N e , which is never zero. N e is defined to be the distance along the ellipsoidal 
normal from the surface of the ellipsoid to the z-axis (see Figure 1). 
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z 



FIGURE 1. ELLIPSOIDAL NORMAL 


From Figure 1, 



Also, 


3. Z 

tan 6 = —r-, from Equation (2) 
b r 


N c = rsec(}) = r^/l + tan 2 <|) = rJl + 


a 4 z 2 

bV 



and 


.4 2 


z" i . a'z" L . a -b 2 


1- ^r b 4 


— ajl + - 


z =ajl 


a 2 z 

sin<b = cos(j)tan(J) = —r— 
b N. 


• <]) = sin 


-i a 2 z 
b 2 N. 


which is the desired expression for the geodetic latitude. 
Also from Figure 1, 


hcos<t> = r 0 - r = r 0 - N e cos<f>. 
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hsin<J) = z 0 - z = z 0 —N e sin<J) 


Multiplying the first equation by cos <}) and the second by sin <() and adding, 


h = r 0 cos<|> + z 0 sin<|)-N e cos <|) +—sin <)) . 


.2x , b . 2 


cos 2 <)) + - 7 -sin 2 <t>: 


2 . a 4 z 2 
tan (})= 4 2 

a 2 z 2 

~b 2 (b 2 -z 2 ) 

i b 2 2 x 

1 + ^-tan (|) 

a 

1+ z 2 

_ b —z 

1 + tan 2 <)> 

““ 2 2 
a z 

1 « o / O ' 


b 2 (b’-z 2 ) 


b 4 +(a ! -b 2 )z ! 


Hence 


1 a 

/2 2 = —using Equation (13). 
e z N: 


h = r 0 cos<(> + z 0 sin <(> - — , 


which is the desired expression for the geodetic altitude. 

To recapitulate, given a, b, xo, yo, and zo, the algorithm is as follows 

7 - b ,7 a / t " T l^nl 7 


e 2 _,b_ ,2 _ a_ _ / 2 , 2 N 2 ,2 r p 

6—1 2 ,e t.2 b r 0 V X 0 + y 0 . P ,2 » 9 P D + 2 /2 ‘ 

a. D 6 6 6 


If q > 0, then 


2 2 

p= 27aXgl Q = M P + q 3 + Vpp, Q' = (^/P + q 3 - VF) 1 , t = ^(q + Q + Q'), 


t, z==—^c-~p + *^2p 2 - q-2t-2pc + 4^/t^Tp^b^Jsign z 0 , 


e' 2 z 2 _ a 2 z a 2 

N e = a. 1 + — 2 ~, <t> = sin 1 ~ 2 —, h = r 0 cos<t> + z 0 sin<t>-— . 
V b bN p N„ 
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CONCLUSION 


An efficient model for transforming Earth-centered, Earth-fixed coordinates to geodetic 
coordinates has been presented. The model systematically derives exact expressions for the 
geodetic latitude and altitude, which are free of singularities. To the author’s best knowledge, no 
such expressions exist in the literature. 


10 


NSWCDD/TR-94/77 


DISTRIBUTION 

COPIES COPIES 

INTERNAL 


DOD ACTIVITIES (CONUS) 

ATTN SP231 1 

STRATEGIC SYSTEMS PROGRAMS 
DEPARTMENT OF THE NAVY 
WASHINGTON DC 20376-5002 

DEFENSE TECHNICAL INFORMATION 
CENTER 

CAMERON STATION 

ALEXANDRIA VA 22304-6145 12 

ATTN CODE E29L (TECHNICAL LIBRARY) 1 

COMMANDING OFFICER 

CSSDD NSWC 

6703 W HIGHWAY 98 

PANAMA CITY EL 32407-7001 


NON-DOD ACTIVITES (CONUS) 

ATTN GIFT AND EXCHANGE DIVISION 4 

LIBRARY OF CONGRESS 
WASHINGTON DC 20540 

THE CNA CORPORATION 
PO BOX 16268 

ALEXANDRIA VA 22302-0268 1 

ATTN RICHARD RAPP 1 

DEPARTMENT OF GEODETIC SCIENCE 
AND SURVEYING 
404 COCKINS HALL 
1958 NEIL AVE 
COLUMBUS OH 43210-1247 


E231 


3 

E282 

SWANSBURG 

1 

K10 


3 

K104 

FELL 

1 

K12 

CARR 

1 

K12 

O’TOOLE 

1 

K13 

EVANS 

1 

K407 

GATES 

1 

K41 

SOFAIR 

10 

K42 


1 

K43 


1 

K44 


1 

K52 


1 

N74 

(GIDEP) 

1 


(1) 



REPORT DOCUMENTATION PAGE 


Form Approved 
OBM No. 07040188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time fa reviewing instructions, search existing data sources, gathering and 
maintaining the data needed, and completing and reviewing the collection of information. Send oomments regarding this burden or any other aspect of this collection of information, including 
suggestions for reducing this burden, to Washington Headquarters Services. Directorate for information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 

22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington, DC 20503. 

1. AGENCY USE ONLY ( Leave blank) 

2. REPORT DATE 

April 1995 

3. REPORT TYPE AND DATES COVERED 

Final 

4. TITLE AND SUBTITLE 

Improved Method for Calculating Exact Geodetic Latitude and Altitude 

5. FUNDING NUMBERS 

6. AUTHOR(s) 

Isaac Sofair 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Commander 

Naval Surface Warfare Center, Dahlgren Division (Code K41) 

17320 Dahlgren Road 

Dahlgren, VA 22448-5100 

8. PERFORMING ORGANIZATION 

REPORT NUMBER 

NSWCDD/TR-94/77 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

10. SPONSORING/MONITORING AGENCY 

REPORT NUMBER 

11. SUPPLEMENTARY NOTES 

12a. DISTRIBUTION/AVAILABILITY STATEMENT 

Approved for public release; distribution is unlimited. 

12b. DISTRIBUTION CODE 

13. ABSTRACT (Maximum 200 words) 

An efficient model for transforming Earth-centered, Earth-fixed coordinates to geodetic coordinates is presented. The model 
systematically derives exact expressions for the geodetic latitude and altitude, which are free of singularities. To the 
author’s best knowledge, no such expressions exist in the literature. 

14. SUBJECT TERMS 

Geodetic Latitude, Geodetic Altitude 

15. NUMBER OF PAGES 

18 

16. PRICE CODE 

17. SECURITY CLASSIFICATION 
OF REPORT 

UNCLASSIFIED 

18. SECURITY CLASSIFICATION 
OF THIS PAGE 

UNCLASSIFIED 

19. SECURITY CLASSIFICATION 
OF ABSTRACT 

UNCLASSIFIED 

20. LIMITATION OF ABSTRACT 

UL 


NSN 7540-01 -280-5500 Standard Form 298 (Rev 2-89) 


Prescribed by ANSI std. Z39-18 
298-102 






GENERAL INSTRUCTIONS FOR COMPLETING SF 298 

The Report Documentation Page (RDP) is used in announcing and cataloging reports. It is important that this 
information be consistent with the rest of the report, particularly the cover and its title page. Instructions for filling 
in each block of the form follow. It is important to stay within the lines to meet optical scanning 
requirements 


Block 1. 


iiencv Use Only (Leave blank). 


Block 2. Report Date . Full publication date including 
day, month, and year, if available (e.g. 1 Jan 88). Must 
cite at least the year. 

Block 3. Type of Report and Dates Covered . State 
whether report is interim, final, etc. *lf applicable, enter 
inclusive report dates (e.g. 10 Jun 87 - 30 Jun 88). 

Block 4. Title and Subtitle . A title is taken from the 
part of the report that provides the most meaningful 
and complete information. When a report is prepared 
in more than one volume, repeat the primary title, add 
volume number, and include subtitle for the specific 
volume. On classified documents enter the title 
classification in parentheses. 

Block 5. Funding Numbers . To include contract and 
grant numbers; may include program element 
number(s), project number(s), task number(s), and 
work unit number(s). Use the following labels: 


C - Contract 
G - Grant 
PE - Program 
Element 


- Project 

- Task 

- Work Unit 
Accession No. 


Block 6. Author(s) . Name(s) of person(s) responsible 
for writing the report, performing the research, or 
credited with the content of the report. If editor or 
compiler, this should follow the name(s). 

Block 7. Performing Organization Name(s) and 
address(es) . Self-explanatory. 

Block 8. Performing Organization Report Number . 

Enter the unique alphanumeric report number(s) 
assigned by the organization performing the report. 

Block 9. Sponsoring/Monitoring Agency Namefs) and 
Address(es) . Self-explanatory. 

Block 10. Sponsoring/Monitoring Agency Report 
Number . (IfKnown) 

Block 11. Supplementary Notes . Enter information not 
included elsewhere such as: Prepared in cooperation 
with...; Trans, of To be published in... . When a 
report is revised, include a statement whether the new 
report supersedes or supplements the older report. 


Block 12a. Distribution/Availability Statement . 
Denotes public availability or limitations. Cite any 
availability to the public. Enter additional limitations 
or special markings in all capitals (e.g. NOFORN, REL, 
ITAR). 


DOD - See DoDD 5230.24, “Distribution 

Statements on Technical Documents” 
DOE - See authorities. 

NASA - See Handbook NHB 2200.2 
NTIS - Leave blank 


Block 12b. Distribution Code. 


NASA 

NTIS 


Leave blank. 

Enter DOE distribution categories from 
the Standard Distribution for 
Unclassified Scientific and Technical 
Reports. 

Leave blank. 

Leave blank. 


Block 13. Abstract . Include a brief (Maximum 200 
words) factual summary of the most significant 
information contained in the report. 

Block 14. Subject Terms . Keywords or phrases 
identifying major subjects in the report. 

Block 15. Number of Pages . Enter the total number 
of pages. 

Block 16. Price Code . Enter appropriate price code 
(NTIS only). 

Block 17-19. Security Classifications . Self- 
explanatory. Enter U.S. Security Classification in 
accordance with U.S. Security Regulations (i.e., 
UNCLASSIFIED). If form contains classified 
information, stamp classification on the top and 
bottom of this page. 

Block 20. Limitation of Abstract . This block must be 
completed to assign a limitation to the abstract. 

Enter either UL (unlimited or SAR (same as report). 

An entry in this block is necessary if the abstract is to 
limited. If blank, the abstract is assumed to be 
unlimited. 


Standard Form 298 Back (Rev 2-89) 





